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Abstract. We study the ratchet effect of a damped relativistic particle driven by 
both asymmetric temporal bi- harmonic and time-periodic piecewise constant forces. 
This system can be formally solved for any external force, providing the ratchet 
velocity as a non-linear functional of the driving force. This allows us to explicitly 
illustrate the functional Taylor expansion formalism recently proposed for this kind of 
systems. The Taylor expansion reveals particularly useful to obtain the shape of the 
current when the force is periodic, piecewise constant. We also illustrate the somewhat 
counterintuitive effect that introducing damping may induce a ratchet effect. When the 
force is symmetric under time-reversal and the system is undamped, under symmetry 
principles no ratchet effect is possible. In this situation increasing damping generates 
a ratchet current which, upon increasing the damping coefficient eventually reaches a 
maximum and decreases toward zero. We argue that this effect is not specific of this 
example and should appear in any ratchet system with tunable damping driven by a 
time-reversible external force. 
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1. Introduction 

The ratchet effect is identified with the motion of particles or solitons induced by 
zero-average periodic forces jH sometimes in the presence of thermal fluctuations. 
The effect arises as a subtle interplay between nonlinearities in the system and broken 
symmetries. Ratchets appear in many fields of physics, where net motion is generated 
either by an asymmetric, periodic, spatial potential [3j SI El d El HD], or by an 
asymmetric temporal forcing [TQl EH EEJ E31 dU [El EH CEO HE] • In both cases the ratchet 
effect can be regarded as an application of Curie's symmetry principle, which states that 
a symmetry transformation of the cause (forces) is also a symmetry transformation of 
the effect (ratchet velocity) [19} 120]. 

Most studies of ratchets driven by temporal forces employ a bi-harmonic forcing 



where p and q are positive integers which, without loss of generality, can be taken co- 
prime (otherwise common factors can be absorbed in the frequency u) and e±, €2 are 
small non-zero parameters. If both p and q are odd, the force (Cp) exhibits the shift 
symmetry (Sf)(t) = f(t + T/2) = —f(t), where T = 2tt/u. In systems invariant under 
time translations this implies that both, f(t) and —f(t) generate the same ratchet 
current (or velocity) defined as [211 122] 



where x(t) is the position of the particle, soliton, or localized structure. If reversing 
the force changes the sign of the current, this current must be zero. So shift-symmetric 
bi-harmonic forces cannot induce a ratchet effect. In contrast, if p and q have different 
parity, shift symmetry is broken by f(t) so the force can induce a nonzero net current 



Many attempts have been made to determine quantitatively the dependence of the 
ratchet velocity, v, on the parameters of the bi-harmonic force (0Q) (TTJ EH |25] . Invariably, 
the analysis performed in these works rests on the so-called method of moments, where 
it is assumed that the average ratchet velocity can be expanded as a series of the odd 
moments of f{t), i.e. EfcLi(Lf(*)] 2fc+1 ) with (M*)) = Iq dth{t). This method seemed 
to work for some systems but not for others without a clear reason and with no known 
criterion to tell ones from the others. We have recently shown that the moment method 
relies on an assumption that almost never holds, and have provided an alternative 
procedure that yields the correct result regardless of the system [23] . 

The aim of this paper is to provide explicit examples which illustrates this otherwise 
abstract method — the functional expansion of v in terms of / — using a working example 
for which an analytic solution can be found. The system represents the motion of a 
damped, relativistic particle under the effect of two different forces: a bi-harmonic force 




(1) 



v= lim - / x(r)dr= lim 





(2) 



[23J 
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like ([T]), and a time-periodic piecewise constant force like 

r e x if o < t < t u 

/(*)=< if 7} <t<T-T,, (3) 

[ -ei if T — Ti <t <T. 

To this purpose we introduce the model as well as its analytic solution in section [2j 
In section [3] we discuss the phenomenon of damping-induced ratchets. The formalism 
developed in [23] is fully illustrated for this problem in section |U For these two specific 
driving forces it is also shown that the method of moments is valid only when the 
dynamics of the relativistic particle is overdamped, and fails otherwise. Conclusions are 
summarized in section |5j 



2. Motion of a relativistic particle driven by a bi-harmonic force 

The equation of motion of a relativistic particle with mass M > 0, whose position and 
velocity at time t are denoted x(t) and u(t), respectively, is 

37 = u(t), x(0) = x , (4a) 

at 

M Tt = ~ fm ~ u " )3/2 ~ 7M(1 " m2) ' m(0) = Mo ' m 

where xq and u$ are the initial conditions, 7 > represents the damping coefficient 
and f{t) is a T-periodic driving force. Notice that if the force fit) satisfied (Sf)(t) = 
f{t+T/2) = —fit), then ( jjjjj) would be invariant under a combination of shift symmetry 
(<S : 1 1— > t + T/2) and the sign change x —x. 
Changing the variable u{t) by the momentum 
= Mu( t) 

- «m 

transforms (jjfej) into the linear equation 

§ = -l,P- m P(0) = P = ^£i=, (6) 
at y/i - u 2 

where (3 = 7/M. Equation (JSJ) is easily solved to give 

Pit) = Poe'^ - [ dzfiz)e-^- z \ (7) 
J 

From (jSJ) one obtains 

^flV (2k-l)\\ fPit)\ 2k+1 

Let us now focus our attention on the T-periodic driving force f{t) given by (CQ) with 
p = 2 and q = 1 (the most common choice of parameters E21 Ell I2"7|). Substituting 
dU into © leads to 

P(t) = -Poe _/3< - ei cos(wt + 0i - Xi) - h cos(2o;i: + 2 - X2), 
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with 

P = P + ?i cos(0i - xi) + e 2 cos(0 2 - X2), 
h =e 1 (p 2 + co 2 )- 1/2 , Xi = tan- 1 (u/fi) , 

e 2 = e 2 (/3 2 + 4a; 2 )" 1 / 2 , X 2 = tan^ 1 (2w/0) . 

As t — > 00 the momentum P(t) behaves, for any /3 > 0, as 

P(t) ~ — ei cos(w£ + 0i — Xi) - e 2 cos(2u;t + 2 - X2), 

thus the term P(t) 2k+1 in (J8j) is 0(e\e?f) with r + s = 2k + 1. Since the time average of 
P(t) is zero, the leading term of in powers of €\ and e 2 will be 

~2lP \ I P(T)3 dT = 2W e1e ~ 2 / C ° S(WT + 01 " Xl)2 COs(2wr + 02 " X2) dT 

3 

r e 2 e 2 cos(20i - 2 + % 2 - 2%i). 



8M 3 

Therefore, the rachet velocity (jSJ), for small amplitudes ei and e 2 , is given by 

v = Bt\t 2 cos(20! - 2 + O ) , (9) 

with 

3 / 2uj 3 \ 

B = ; = X2 - 2yi = - tan" 1 — — — , (10) 

8M 3 (/3 2 + o; 2 )#+4? A A V/3(/3 2 + 3w 2 ); 1 ; 



in agreement with the result reported in 

Notice that in the undamped limit 7 — > (equivalently — >■ 0) the parameters ffTUl) 
become 

R = 3 fl = JL 

16AW 2' 

whereas in the overdamped limit M — > (and therefore /3 — ?• 00 with M/3 = 7) we get 

B = ^ ^ = ' 
both limits agree with the predictions of [23J. 

3. Ratchet induced by damping 

The depence of v on parameters of the system like the damping coefficient (through 
= 7/M) shown in ([9]) and (flOl) reveals an interesting effect. If we take <f\ = </> 2 = in 
f(t) and do some algebra, the ratchet velocity (for small amplitudes of the force) turns 
out to be 

^VtfM V{x) = , ^fjf, ^ . (11) 



oo 3 M 3 ^' " v ' 8(x 2 + l) 2 (x 2 + 4)' 

Function V(x) is depicted in Figure [TJ The most remarkable observation is that the 
current increases up to a maximum with increasing damping before it begins to show 
the expected decay. Intuition dictates that the current should decrease with damping 
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Figure 1. Plot of the current velocity v, in units of e\e2 / '{M w) 3 , vs. the damping 
coefficient j3, in units of the frequency, u>, induced by a biharmonic force like ([1} with 
<t>i = 4>2 = 0. Notice that this force is time reversible, i.e., f(—t) = f(t). 



because friction opposes movement, so the fast increase it reveals for small damping is 
counterintuitive. 

The cause of this effect is the interplay between the breaking of the time-reversal 
symmetry 1Z : t y — t that generates the ratchet current, and the damping that hinders 
it [27]. In the limit /3 — > the system (j2J) is invariant under 71 and a sign change of 
u, because for 4>i = 02 = the force (QQ) satisfies f(—t) = f(t). Accordingly v = in 
this limit. But introducing damping breaks the symmetry of the equation and induces 
a net movement of the particle. For small damping, the higher the damping coefficient 
/3 the larger v. If we keep on increasing /3 eventually the friction it introduces in the 
movement of the particle causes the decay of v as /3 -3 . 

This argument makes it clear that in any ratchet system with a tunable damping 
and undergoing the action of a time-reversible bi-harmonic force, the ratchet effect can 
be generated upon increasing damping above zero. 

4. Ratchet velocity as a functional of the force 

The starting point to obtain formula (Q for a ratchet system is to realize that v is a 
functional of f(t) and that, under certain regularity assumptions, one such functional 
can be expanded as a functional Taylor series [281 1291 130] as 



where the kernels c n (ti, . . . , t n ) are proportional to the nth functional derivatives of 
the functional v[f}. These kernels can be taken T-periodic in each variable and totally 
symmetric under any exchange of variables. Only odd terms appear in this expansion 
as a cosequence of the symmetry v[—f] = —v[f] that these systems have. 




(12) 
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That v is indeed a functional of f(t) in this example is obvious from equations ([2])- 
([7]). The aim of this section is to determine explicitly the expansion f lT2|) for this exactly 
solvable example. 

Let us start off by rewriting the integral in ([7]) as 

t dzf(z)e-^ = h(t) + I 2 (t), (13) 

n(t) 

-p{t-z-(k-l)T) 



fe=l ^ 
Mi) 

I 2 (t) = / dz/(^)e- /3(a( * ) - z) ; 



where a(t) = t — n(t)T and n(t) = [t/T] ([X] denoting the integer part of X). Notice 
that a(t + T) = a(t). Now, since 

s(t) ^^^-^ = i^-i, (14) 

fc=l 

then h(t) = e- pt CS(t), with 

C= f 'dzf(z)[e/* - 1]. (15) 
./o 

Using this form in (J7J) we can write 

P(t) = Ae- /3 * + P(t), (16) 
where A = P(0) + C(e^ T — 1) _1 and P(t) is the T-periodic function 

P W = — Sf-^T / dy/^K^V - 1] - / dj//(y)e-« a <*>-*>. 
e — 1 Jo Jo 

It is thus enough to obtain P(t) in the interval < t < T, where it adopts the compact 
form 

P(t) = - f dyf(y)e-^ X (y,t), (17) 
Jo 

defining 

1 _ Q-PV 

x(y,t)= ef}T _ 1 +Q(t-y) (18) 

(as it is customary, 6(x) denotes the Heaviside function, which is 1 if x > and 
otherwise). 

Equations (JTTj) - (fT8"]) have a well defined (3 — > + limit, namely 

P{t) = ~f dyf(y) X i{y, t), X i(y, t) = ^ + Q(t-y). (19) 

On the other hand, for zero- average forces f(t) the kernel x{t, z) can be further simplified 
to 

P(t) = ~ f dyf(y)e-^ X 2(y,t), Xa (y,f) = —L^ + 6(f - y). (20) 
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Whatever the form, it should be periodically extended beyond the interval [0,T). 
It is then clear that ([2]) and (jSJ) boil down to 

A direct comparison of this equation with the functional Taylor series (|T2j) yields 

P2fc(ti,...,t 2 fc) = 0, (22a) 

C2fc+i(ti, • • • ,t 2 k+i) — (— ^ J -j^-j^^T 2k a 2 k+i(ti, ■ ■ ■ ,t 2 k+i), (226) 

where 

a m (t 1 ,...,t m )= dre-^ r -*")TTx(t fe ,r), t=-Vt (23) 

J ° L\ m tl 

As expected [23], functions a m (ti, . . . , t m ) are, by construction, T-periodic in each 
variable and symmetric under any exchange of their arguments. 

The integral in (|23|) can be performed integrating by parts and taking into account 
that -^x{Vi T ) = ^( r ~ v) ( a Dixac delta). The result is 



in 



a m (t 1 ,...,t m )=—in x (t k , o)-n x(t k +t,o)+yi ^ mt] n ^ r > ( 24 ) 

^ lfc=i fc=i i=i fc=i, fe^j J 

where we have used the fact that xi^k, T)e _/3T = x(tk + T, 0). As usual, empty products 
are assumed to be 1 (the case of the last term for m — 1). 

The limit /3 — ?• of this expression is better obtained by replacing x(y, t) by XiiVi t) 
in ( 1231) and integrating by parts again. This results in 



a m (ti, . . . ,t m ) = T Y[xi(tk,T) - ^tj Yl Xi(tk,tj). (25) 

fc=l j=l k=l, te£j 

Finally, in the overdamped case (M — > 0, /3 — > oo), instead of the evolution of 
P is given by P(t) = — (l//3)/(t), so v can be expressed simply as 

From (fT2l) and (126]) it follows that C2fc(ti, ■ ■ ■ ,t 2 k) =0 and 

/ y2\ ^fc - 1)!! 

c 2 fc+i(*i, • ■ -,t 2k+ i) = - ( — —J 2A;+1 " ffii - fr») • • • - t 2fc+ i). (27) 
^.i. Forcing with a time-periodic piecewise constant force 



The expansion (fT2l) with kernels (122 £>j) and ( 123]) turns out to be useful to analyze different 



types of forcing. For instance, another standard choice in the literature (see pQ and 
references therein), alongside with the bi-harmonic force (CQ), has been the time-periodic 
piecewise constant force defined in ([3]). This force is shift-symmetric only for T} = T/2, 
so any other value T/ < T/2 breaks this symmetry and induces a ratchet current. 
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In order to ascertain the effect of this force in system (|46|) for small amplitudes 
ei <C 1, we will compute the first nonzero term in the expansion ffl2|) . To that purpose 
we need to evaluate (c.f. equation ( 1231 ) 

K m = (a m (£ 1; . . . , O/(*0 • • • /(O) = / T [e^ T /(r)] m dr, (28a) 

Jo 

J(t)ee1 / e^ X2 (t,r)f(t)dt, (28b) 



T „ 

where the choice X2(^, r ) instead of r ) is made because fit) in d3J) has zero average. 
According to (1201 X2(^?t) = (1 — e~ /3T )~ 1 X2(£j t), where 

,1 if t < r, 

&(*.t)H -/3T _ (29) 



e- pi if t > 



r. 



Substitution into (128 fej) yields 



/(r) = ^ 



1 sinh 2 ^ 1-0(7-) 



e 



-/w V 2 



(30 a) 



e pr _ e /3T ; if < t < Tj, 

{ if/, < v < V 7,. (306) 

e «T-T,)_ e ^r if T-Tt<T<T. 

It is straightforward to check that l£i in (128 a\i vanishes, so the first term that may 
not be zero is K3. Lengthy calculations lead to 



Ki = ~w w^w sinh I — 2 — J sinh It" J ' (31) 

that is to say 

- ^ ^( p{T : 2T ' ) Y n K im/2) ^o(4). (32) 



((3Mff3T V 2 ; sinh 2 (/3T/2) 

It is interesting to noticing that K3 = if T\ = T/2 because in that case the time- 
periodic piecewise constant force ([3]) is shift-symmetric. On the other hand, we can 
determine the value of T) for which the ratchet effect is maximum by differetiating (|32|) . 
This leads to 

sinh (^- 3r <>) smh smh 3 = o. (33 ) 

The only three solutions to this equation are 7] = 0, 7} = T/2 and 7} = T/3. The first 
two do not produce any ratchet current (with 7] = / = whereas for 7} = T/2 the 
force is shift-symmetric), therefore the last one provides its maximum value. 

As a final remark, expression ( l32l has well defined overdamped (M — > 0, /3 — > 00, 
with finite 7 = (3M) and undamped (/3 — > 0) limits. In fact, the undamped limit of (1321) 
yields 

T 4 (T _2T / ) 2 + o(e?), (34) 



4(MT) 
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whereas the overdamped produces v = o(ef). Indeed, since f(t) 2k+1 = el k f(t), in the 
overdamped case Eq. (1261) immediately implies v — 0. This is in marked contrast with 
the overdamped deterministic dynamic of a particle in a sinusoidal potential driven by 
a bi-harmonic force [31]. In this case, the zero ratchet velocity can be explained as a 
symmetry effect. Indeed, notice that fit) = —f(—t) when f{t) is given by (131) (something 
that only happens for the bi-harmonic force ([1]) for specific choices of the phases), and 
that the overdamped limit of Eq. (|4fej) remains invariant under a simultaneously action 
of time-reversal and a sign change of x and u (see [23J for further details). 

5. Discussion 

We have studied the dynamics of a damped relativistic particle under two zero 
average T-periodic forces which breaks the shift-symmetry fit + T/2) = —f{t). This 
nonlinear system can be explicitly solved through a transformation that renders it linear. 
Therefore, the ratchet average velocity, v, is exactly obtained for any arbitrary force f{t). 
This result allows us to show, first of all, that the ratchet velocity cannot be obtained in 
general by using the method of moments (according to which v is obtained as a series of 
the odd moments of fit)). And secondly, that v is a functional of f(t), i.e. v[f}. Indeed, 
for any T-periodic force we have explicitly found the coefficients of the functional Taylor 
expansion ([121) . In particular, this expansion shows that the method of moments is only 
justified in the strict overdamped limit (see Eqs. ( 1261) and (1271) ). Due to the symmetry 
v[f] = —v[—f] only odd terms contribute to the Taylor expansion. Besides, since the 
ratchet velocity is translationally invariant, the kernel ci(ti) must be a constant. So the 
first order term vanishes because the force has zero average. Therefore, the first term in 
the expansion (fl2l) that is not necessarily zero is the third one, irrespective of the kind 
of nonlinearity of the system. 

We have chosen to illustrate this functional representation the bi-harmonic force ([1]) 
(with p = 2 and q = 1) as well as a time-periodic piecewise constant force ([3]). We have 
obtained the leading term of the average velocity for both these forces. They are given 
by equations (l9l)- (TIOl) and ( |32l) . respectively. It is worth emphasizing that the method 
of moments always predicts a zero ratchet velocity when the system is driven by a time- 
periodic piecewise constant. This is to be compared with the result f[32]) obtained here. 
We have discussed the two limiting dynamics: undamped and overdamped. In these two 
limits the system remains invariant if the driving force has the symmetries f(t) = f(—t) 
and f(t) = — /(— i), respectively. If the relativistic particle is driven by a bi-harmonic 
force v ~ efe 2 cos(20! — </> 2 ) in the overdamped limit, whereas v ~ e\e 2 sin(2</>i — 2 ) in the 
undamped limit. In the latter case this means no ratchet current if we set 0i = 2 = 0. 
The unexpected consequence of this is that introducing damping generates a ratchet 
current, whose intensity grows up to a maximum before it drops to zero upon a further 
increase of the damping. This effect is a result of a trade off between symmetry effects 
and friction and our prediction is that should be observed in any system with damping 
and forced with a time-reversible external force. 
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On its side, if an overdamped relativistic particle is driven by a time-periodic 
piecewise constant force like (J3J), the ratchet velocity is always zero as a consequence of 
the symmetry f(t) = —f{—t) exhibited by this force. 

Summarizing, we hope to have illustrated the predictive power of the Taylor 
functional expansion method introduced in [23J. This working example also shows 
that this is the only reliable method to analyze the ratchet current as a function of 
the parameters of the external force. The most widely used alternative so far, the 
method of moments, is shown to work only in the overdamped limit of the dynamics of 
a relativistic particle driven by a periodic force. When damping is finite and the forcing 
of the system is bi-harmonic, the rathet current predicted by the method of moments 
still retains some relevant features of the exact one (Q. However, it dramatically fails 
if the system is driven by a piecewise constant force, because it always predicts a zero 
ratchet current, in marked constrast with the result ( |32l) predicted by the functional 
Taylor expansion. 
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